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X. INTRODUCTION 


The progress achieved In the third year of this three-year research 
program is reported. The overall research program is organized into three 
different subjects: (1) the Jupiter system, (2) the Saturn system, and 

(3) Comets. In the Jupiter system, the extended atmospheres of Io and their 
relations and interactions with the planetary magnetosphere are the topic of 
interest. In the Saturn system, the hydrogen torus of Titan, the atmospheres 
associated with the planetary rings, and the extended atmospheres of the 15- 
ring icy satellites as well as the planet are the topics of concern. For 
comets, the dust and gas clouds in the inner and outer coma are the areas of 
focus . 

The three primary goals of the overall program are (1) to better under- 
stand the characteristics of the local atmosphere of Io and the manner in 
which it interacts and delivers a significant source of heavy-ion plasma to 
the planetary magnetosphere, (2) to understand the spatial distribution of 
hydrogen in the Saturn system in terms of the planet, its rings and its 
satellites, and to better define the source of protons for the magnetosphere, 
and (3) to provide fresh insight into the chemical and physical nature of 
comets through study of their extended atmospheres. These research objectives 
are important in achieving better understanding of the basic satellite and 
magnet os phuric processes in the outer solar system and basic comet and solar 
radiation processes in the inner solar system. To pursue these goals, the 
research program has employed both exploratory calculations (to uncover impor- 
tant physical mechanisms) and more-refined model calculations to evaluate the 
quantitative significance of these mechanisms. 

For the Jupiter system, the overall three-year research plan is summa- 
rized in Table 1. There are three major subjects. The first subject, the Io 
sodium cloud, is divided into five subtopics. Efforts during the third year 
have been focused upon subtopics (i), (iii) and (iv). Efforts originally to 
be applied to subtopic (ii) were reassigned to topic (iv) where the same 
objective was able to be pursued more effectively using the more recently 
discovered directional feature data of Pilcher (1982). For the Io potassium 
cloud, the second topic in Table 1, model development was completed this year 
and applied to observational data of Trafton (1981). The third topic of Table 
1 has remained inactive, since no new gas clouds were discovered during the 
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JUPITER SYSTEM: TIIREE-YEAR PLAN FOR MODEL DEVELOPMENT AND ANALYSIS 
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third year. Progress achieved in the first two subjects of Table 1 is 
discussed in section II. 

For the Saturn system, the overall three year research plan is summarized 
in Table 2 and is organized into the study of font objects, each known or 
expected to provide a source of hydrogen atoms and ultimately protons to the 
circumplanetary space. During the third year, major emphasis has been placed 
upon better understanding the role of Titan (and lienee ultimately non-Titan 
sources) in the new light that the planet, may itself also provide a signifi- 
cant source of hydrogen to the Saturn system (Shomanskv and Smith, 1982). 
Progress resulting from these efforts is discussed in Section III. 

For comets, our research efforts have been reactivated in the third year 
as scheduled in Table 2. The AER particle trajectory model has been further 
developed and applied to the extended hydrogen atmospheres of comets. Results 
of these, efforts are reported in section IV. 
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It. TUB .HWITBR SYSTEM 


2 . 1 In t r> Hjiu't iim 

The th*vo-y<*,u research pi. inn for the To sodium eJoml and Io potassium 
el oml ar»» summarized in Table 1# ‘'todies tor the Io sodium eloud are subdi- 
vided into a general model development and four separate subtopics organized 
around the different types ot observational data that are available. Progress 
in the third year has been built upon our primary model-development activities 
ot the lirst year and the initial phase ot model application undertaken during 
tile se. ond year. The improved models, which now include the oscillating Io 
plasma torus as a time-space varying sink lor sodium and potassium as well as 
evioi more complex space-time dependences, provide new and powerful modeling 
tools that have al lowed us to draw some significant conclusions from observa- 
tional data. These conclusions are discussed below. 

2 . 2 I n t e r a et:lj)ns _<)f the I o S < >< 1 lie n Cloml and PI nsma _ Torns^ 

The space-time history of an ensemble of sodium cloud atoms initially 
ejected from Io and moving through the inner magnetosphere of Jupiter depends 
upon the physical processes by which these atoms may interact with the plasma 
torus. There are four basic processes: electron Impact ionization, charge 

exchange reactions, photoionization, and elastic ion-neutral collisions. As 
discussed below, ionization by electrons lppears to be by far the dominant 
process for shaping the spatial morphology of the sodium cloud in Region B and 
in near Region G. 

The first three of the plasma-atom interactions noted above effectively 
remove sodium from the cloud and hence determine its lifetime in the plasma 
torus. Ionization of sodium by electron impact ionization Is very rapid in 
the plasma torus ( 1 hour lifetime based upon Voyager l plasma data; see 

Figure 1). Charge exchange reactions remove neutral sodium from the cloud and 
create in its place either an ion (e.g., X ++ + Na + X + + Na + ) or a fast 
neutral (~60 km sec -i , which was the corotating ion X + before the reaction 
X + + Na + X + + Na + ). The fast neutral so produced rapidly escapes from the 
Jupiter system unless it if: trapped in the outer magnetosphere. For the case 
of resonance charge, exchange (X^ = Na + ), the fust neutral is sodium, hut it 
moves so rapidly that it appears to contribute minimally to the bright Region 
B and close Region C portion of the cloud, although it may be a very important 
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VERTICAL DISPLACEMENT FROM THE CENTRIFUGAL EQUATOR (Rj) 



Figure 1. Sodium Electron Impact Ionization Lifetime in the lo Plasma Torus. 

The lifetime calculation from Smyth (1983) is based upon the. cross- 
section measured by McFarland (1965), McFarland and Kinney (1965), and 
Zapesochnyi. and Aleksnkhin (1968), Voyager 1 in situ measurements of 
the plasma density (Bridge, Sullivan and Bagenal, 1980; Bagenal and 
Sullivan, 1981), and Voyager 1 UVS (Shemansky, 1980) and JLn_ situ 
measurements (Scudder, Sittler and Bridge, 1981) for the electron 
temperature. The two-dimensional common temperature model of Bagenal 
and Sullivan (1981) was selected to describe the electron density. 
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source *'or the fast component of the fainter remote-sodium observed by Brown 
and Se'mv.iwor (1981) and more recently near lo by Schneider, Hunten aid Brown 
( 19 ^ 4 ) and Tr.iuger (1984) Charge exchange reactions for Region B and close 
Region C may therefore be treated simply as a lifetime process. Based upon 
the best Voyager 1 data (Bagenal 1984; Shematisky, 1984) and available cross- 
sections (Johnson and Strobel, 1982), the charge, exchange lifetime for sodium 
in the plasma torus is dominated, by more than an order of magnitude, by 
electron impact ionisation for radial distances beyond about 5.5 Rj (where tae 
electron temperature „> 4 ev), and hence charge exchange becomes ineffective. 
Inside of about 5.? Rj, the electron arid ion temperatures become cool, the 
electron impact ionization lifetime becomes consequently larger (see Figure 1) 
and the plasma torus collapses latitudinally to within about 0.2 Rj of the 
plasma symmetry plane. The lifetimes for resonance charge exchange and for 
non-resonance charge exchange (only 0 ++ , S ++ and S +i + have significant cross- 
sections with Na) therefore eventually compete and become more important than 
electron impact ionization Inside of 5.5 Rj. At that point, however, the ions 
are confined so closely to the centrifugal symmetry plane of che plasma torus 
and have sufficiently large lifetimes that they are incapable of having any 
significant spatial overlap and effect during the nominal 20 hours of flight 
time of sodium atoms in the forward cloud (see Figure 2). Even less important 
is the photoionization lifetime for sodium atoms near Jupiter. It is about 
400 hours and is only Important as a sodium sink far outside of the plasma 
torus. 

Elastic collisions between plasma torus heavy ions and atoms of the 
sodium cloud may also occur. As noted by Brown, Pilcher and Strobel (1983), 
"hard" ion-sodium atom collisions (resulting from the repulsive core of the 
interatomic potential) cause a simple loss of material from the Io’*connected 
cloud and are, in this way, similar to the resonance charge exchange reaction 
discussed above. "Soft" ion-sodium atom collisions (resulting from the 
attractive force produced by the dipole moment induced by the passing ionic 
charge) will result in much smaller momentum transfer that may sever the 
kinematic connection between the initial escape conditions at Io and the 
morphological structure of the near (Region B, close Region C) cloud. This 
effect of soft collisions, if important, could then greatly complicate the 
simpler model description adopted here. To estimate the importance of these 
soft collisions, Brown, Pilcher and Strobel (1983) used the momentum transfer 
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Figure 2. Interaction of the Sodium Cloud and the Plasma Torus. 

The spatially projected overlap in the satellite plane of the Io 
plasma torus with the solar radiation perturbed cloud shape (solid 
boundary with heavy shading) and with the unperturbed cloud shape 
(heavy dashed boundary) is indicated for the diametrically oppo- 
site satellite phase angles. The cloud shape illustrates the 
envelope of sodium atoms after 20 hours of flight time. 
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cross-section based upon the ion-atom polarization potential.. For the large 
polarizability of sodium, however, this potential provides only an over- 
estimation of the strength of the interaction (Johnson, 1984). This occurs 
because the correct potential for sodium at the values of the impact para- 
meters of interest is actually a combination of a comparable repulsive 
potential and attractive polarization potential which partially cancel and 
reduce the value of the true potential and hence its cross-section. The true 
potential for sodium is not presently available. Adopting this polarization 
cross-section to establish an upper bound, the sore elastic collision time for 
sodium Is evaluated in the plasma torus for Voyager 1 plasma conditions and 
shown in Figure 3. The collision time is more than an order of magnitude 
larger than the electron impact ionization lifetime for sodium beyond about 
5.5 Rj. For the Region B overlap beyond about 5.5 Rj, the average collision 
tine is no smaller than 50 to 100 hours in the oscillating torus. Slocg it 
takes a nominal 5 hours for the fraction of the sodium atoms which survive 
ionization by electrons to exit inside of the plasma torus, this 50 to 100 
hour lifetime implies that 5-10% of these surviving sodium atoms will suffer a 
soft elastic collision. Since this is an overestimation, the real effect 
based upon a true potential for sodium must then be at best only a few 
percent. Inside of 5.5 Rj, the situation is even less promising, because the 
ions are confined so closely to the centrifugal symmetry plane of the plasma 
torus that they are incapable of having any significant overlap and effect 
during the remaining nominal 15 hours of flight time of sodium atoms in the 
forward cloud (see Figure 2). 

2 . 3 The Io Sodium Cloud: Region A 

The east-west intensity asymmetry of Region A, discovered and documented 
by Bergscralh et al. (1975 and 1977), was recently attributed by Smyth (1983) 
to the effects of both solar radiation pressure (for very low velocity escape 
of sodium from Io) and the plasma torus sink (for higher velocity escape of 
sodium from Io). The need, noted by Smyth (1983), to further clarify the role 
played by the Io plasma torus sink, has been initially addressed this past 
year. This has been accomplished by using our new Io sodium cloud model, 
which explicitly includes tne oscillating plasma torus sink, to calculate the 
brightness of Region A as viewed through the rectangular viewing slit adopted 
by Bergstralh et al. (1975, 1977). The effects of the ■ •'dilating plasma 


15 



4 5 6 7 8 9 10 11 


RADIAL DISTANCE FROM JUPITER (Rj) 


Figure 3. Elastic Collision Lifetime of Sodium Atoms in the Io Plasma Toru.:. 

The lifetime (in hours) is based upon the polarization lifetime of Brown, 
Pilcher and Strobel (1.983) and Voyager 1 plasma torus conditions (Bagenal, 
1984; Shemansky, 1984). The smallest collision lifetime of about 20 hours 
occurs on the centrifugal equator at a radial distance from Jupiter of 5.7 
Rj, around which unlabeled lifetime contours of 25 and 30 hours are shown. 
The motion of Io (at 5.9 Rj) resulting from the oscillation of the plasma 
torus about the sateliite orbit plane is also indicated. 
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torus sink, the presence of an east-west electric field in the magnetosphere 
(Ip and Goertz, 1983; Barbosa and Klvelson, 1983), and the impact of the 
offset dipole nature of the planetary magnetic field (Acuna, Behannon, and 
Connerney, 1983) have been preliminarily assessed and are briefly discussed 
here. 

The oscillation of the plasma torus about the satellite plane introduces 
a modulation in the brightness of Region A because, of the varying sodium 
lifetime at Io, but it does not introduce any east-west intensity asymmetry 
(if the plasma (vropecties are independent, of system III magnetic longitude). 

An east-west electric field in the magnetosphere has recently been proposed by 
Ip and Goertz (1983) and Barbosa and Kivelson (1983) to explain the dawn-dusk 
asymmetry in the torus ion emission brightness measured by the Voyager UVS 
instrument (Sandel and Broadfoot, 1982; Shemansky and Sandel, 1982). Their 
proposed electric field essentially shifts the plasma torus to a new center 
that is displaced toward the. east by about 0.18 to 0.20 Jupiter radii. For a 
given radial displacement from Jupiter, this provides an overall relative east 
to west shift of the plasma torus properties by twice this amount and gives 
rise to a significant east-west asymmetry in the sodium sink because of the 
sharp temperature gradient inside of lo's orbit. This asymmetry in the sodium 
sink is graphically illustrated by comparison of electron impact ionization of 
sodium at lo's orbital position in Figure 4 (with no electric field) and 
Figure 5 (with an electric field producing a nominal shift of 0.15 Rj). This 
east-west asymmetry in the sink provides an east-west intensity asymmetry in 
the brightness of Region A that is in agreement with the value observed by 
Bergstralh et al. (1975, 1977). This agreement is shown in Figure 6, where 
model calculations incorporating different values of the shift parameter are 
compared with the observations. Inclusion in the model calculations of a 
0.12 Rj offset of the dipole magnetic field toward 149° System III longitude, 
in addition to the east-west electric field, produces the plasma sink illus- 
trated in Figure 7. This additional offset modulates the east-west intensity 
about the mean value established by the east-west electric field as illus- 
trated in Figure 6. This modulation may account for some of the apparent 
scatter present in the intensity data of Bergstralh et al. (1975, 1977). 
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Figure 6. Ease-West Intensity Asymmetry of the Io Sodium Cloud. 

The east-to-west intensity ratio measured by Bergstralh et al. 
(1975, 1977) is shown by the shaded area. Model calculations for 
different values of the plasma torus shift parameter are shown for 
various assumed physical conditions defined in the legend. 
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2 , 4 The Io^S odj um Cl ourt : Re nton B 

Model calculations for the spatial morphology of the Region B sodium 
cloud have been performed where effects of both the oscillating plasma torus 
sink and the acceleration of cloud atoms by solar radiation pressure in the fi- 
ll nes have been included. Model calculations illustrating these effects for 
lo at both eastern and western elongation are shown projected onto the sky 
plane in Figure 8. These calculations assume radial and isotropic ejection of 
sodium from lo. Comparisons of model-calculated and observed (Murcray, 1978) 
l >2 images of the sodium cloud suggest a satellite source rate of 1-2 x 10^ 
atoms sec“* for Region B. In Figure 8, the more elongated sodium cloud at 
western elongation results directly from the combine d effects of orbit pertur- 
bations by solar radiation pressure and the spatially non-uniform lifetime of 
sodium in the plasma torus. This preferential elongation was anticipated by 
Smyth (1983) and has been observed for the sodium cloud (Goldberg et al., 

1980; Goldberg, 1983). The model results in Figure 8 are displayed as they 
would appear above the satellite plane in Figure 9. 

The preferential elongation of the sodium cloud In the west may be 
enhanced, however, by suitably modifying the ejection conditions of sodium 
from lo. This is illustrated in Figure 10, where sodium is ejected radially 
and uniformly only from the band region (defined in the Figure caption). This 
band ejection region is geometrically very similar to that deduced by Pilcher 
et al. (1984) as the source of high-velocity sodium (~20 km sec"^) which can 
explain the space-time variations of the sodium directional features (Region 
C). The initial directions and speeds of sodium escaping lo that form the 
directional features were noted (Pilcher et al. 1984) to be understood in 
terms of a magnetosphei'ic-wind-drlven escape mechanism as discussed in the 
following section. Future modeling of data describing this east-west orbital 
asymmetry in the sodium cloud will allow a more refined description of the 
ejection conditions to be deduced. This Information will then allow us to 
understand better the gas escape mechanism at lo as well as the properties of 
the local satellite atmosphere. 

2 • 5 The lo Sodium Cloud: Region C 

The directional features of the sodium cloud, discovered by Pilcher (see 
Hartline, 1980), extend from Region B into Region C. In the past year, the 
changing north-south directions of these features have been successfully 
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line of sight 


Figure. 9 Model Calculation for the Io Sodium Cloud. 

The calculations are the same two described in Fig. 8, but here 
show the T >2 emission intensity projected normal to the satellite 
plane. Contour values are, from the outside to inside, 0.1, 0.25, 
0.5, 1,2,3 and 5 kiloRayleighs . 
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understood through our model -data comparisons which have involved an aetlve 
collaborative effort with Pilcher. The extremely interesting results of this 
collaborative effort are presented in the paper of Pilcher et al . (1984), 

This paper is accompanied by a companion paper entitled "Ejection of Atoms and 
Molecules from Io by Plasma-ion Impact" submitted by Hieveka and Johnson 
(1984). 

In short, the directional features can result from a source of high velo- 
city (-20 km sec"*) sodium combined with the oscillating neutral sodium sink 
provided by the plasma torus, flit* phase relations'll p between the features' 
directions and Io's magnetic longitude can be understood if escaping sodium is 
initi. Uy directed at near right angles to Io's orbital motion. The initial 
directions and speeds of sodium atoms escaping To to form the directional 
features can he understood in terms of a maenetospheri c-wlnd-dri von escape 
mechanism. The companion paper of Sieveka and Johnson (1984) addresses the 
subject of ejection of sodium atoms from To and provides support for the 
magnetosphoric-wind-drivon escape mechanism. The estimated satellite source 
for the directional features is -1 x 10* atoms sec and is comparable to or 
slightly smaller than the Region B source deduced above. The similarity of 
these two source rates for very different velocity ranges of the ejected atoms 
provides additional diagnostic characteristics (as noted by Sieveka and 
Johnson, 1984) that characterize the magnet os pherie-wi nd-driven escape 
mechanism. 

2 • 6 T he To Potassium Cloud 

The Io potassium cloud has been observed to radiate in the 7665 A and 
7699 A lines by Trafton (1975) and in the 7699 A line by Trauger, Roesler and 
Munch (1976; Munch, Trauger and Roesler, 1976). The radiation is thought to 
arise from solar resonance scattering. Additional observations of Trafton 
(1977, 1981) have shown that the Io potassium cloud exhibits spatial and 
temporal variations similar to those observed for the sodium cloud. These 
observations indicate that the potassium cloud undergoes periodic fluctuations 
in response to solar radiation pressure (producing the east-west asymmetries) 
and in response to the ionizing influence of Jupiter's plasma torus (producing 
the north-south asymmetries). In addition, temporary directional features of 
potassium from Io have also been observed. Modeling improvements for the 
sodium cloud model discussed in the preceding subsection may therefore be 
incorporated directly in the Io potassium cloud model. 
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An To pot. ns, slum cloud model has boon dovolopod this year. It includes 
tho electron impact ionisation lifetime of potassium in the oscillating plasma 
torus (see Figure 11), tho solar spectrum for each potassium line (see Figures 
12 and 13), and tho acceleration experienced by potassium atoms as they under- 
go solar resol nance scattering. The electron impact ionization of potassium 
in Figure 11 is comparable to the sodium lifetime in Figure 1 in both magni- 
tude and spatial character. Both of the potassium lines appear in solar 
absorption features, and the line at 7665 A is almost completely obscured to a 
ground-based observer because of an optically ttuc.K oxygon absorption feature 
present in the earth's atmosphere. To properly model the intensity morphology 
of the cloud or the line profile shape in the 7665 A emission, th». exact shape 
of the solar absorption feature above the earth's atmosphere must be known as 
well as the relative motion between the cloud atoms and the Earth. The compu- 
tation of the unobstructed solar abosrption feature and the earth's atmospher- 
ic transmission function was provided in the first program year by Kuru cz 
(1982), with whom a collaborative effort for this purpose was established. 

The results of his computations are shown in Figures 12 and Figure 13 and have 
been adopted in the model, calculation to determine both the excitation rates 
and the atom accelerations produced by solar resonance scattering. Because of 
the differences in the shapes of the potassium and sodium absorption features, 
the maximum magnitude of the acceleration produced by solar radiation pressure 
is in fact comparable. In modeling, therefore, the Io potassium cloud should 
exhibit east-west asymmetries of similar magnitude to those calculated earlier 
for the sodium cloud. 

Model calculations for the brightness of the Io potassium cloud in the 
7665 A emission line at both eastern and western elongations are shown in 
Figure 14. These calculations include the oscillating plasma torus and assume 
isotropic ejection of potassium at 2.6 km sec - * from Io's exobase. Adopting a 
brightness of about 100±50 Rayleighs near Io as reported by Trafton (1981), a 
potassium flux of 1.2 x 10^ atoms cnT^sec - * (relative to Io's surface) or a 
total source rate of about 5 x 10^ atoms sec - *' is required by the model. The 
results in Figure 14 have used this source rate and are for Io at a system III 
magnetic longitude of 200°. Solar radiation pressure is responsible for the 
shorter cloud in the east and the longer cloud in the west. This value of the 
potassium source rate is about a factor of 20-40 smaller than the sodium 
source rate deduced for Region B and is similar to that expected based upon 
the relative cosmic abundance of the two elements. 
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RIFUGAL EQUATO 


POTASSIUM ELECTRON IMPACT IONIZATION LIFETIME IN THE 

10 PLASMA TORUS 

(Voyager 1 Plasma Data) 



RADIAL DISPLACEMENT FROM JUPITER ( Rj ) 

Figure 11. Potassium Electron Impact Ionization in the Io Plasma Torus. 

The lifetime calculation is based upon the cross-section measured by McFarland 
(1965), McFarland and Kinney (1965), and Zapesochnyi and Aleksakhin (1968), 
Voyager 1 _in_ situ measurements of the plasma density (Bridge, Sullivan and 
Bagenal, 1980; Bagcnal and Sullivan, 1981), and Voyager 1 UVS (Shemansky, 
1980) and _in_ situ measurements (Scudder, Sittler and Bridge, 1981) for the 
electron temperature. The two-dimensional common temperature model of Bagenal 
and Sullivan (1981) was selected to describe the electron density. 
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III. THE SATURN SYSTEM 

3 . 1 Introduct ion 

For the Saturn system, a condensed summary of our three year research 
program is given in Table 2. The four subtopics listed in Table 2 represent 
four plausible sources of hydrogen for the Saturn system. Efforts in the 
first year of this research program were directed primarily to the planetary 
ring atmosphere, where new data had become available from the Voyager 1 and 
Voyager 2 encounters. In the second year, the planet was recognised as being 
a possible source of K for the larger Saturn system (Shemansky and Smith, 

1982), and the scope of the problem has been appropriately broadened and 
refocused as indicated by the addition of the first topic in Table 2. An 
extended hydrogen atmosphere for Saturn must now be considered as a plausible 
source of H atoms previously interpreted as a planetary ring atmosphere and 
might even be an important source for some of the hydrogen previously asso- 
ciated with a Titan source. 

To further explore and understand this subject, our past collaborative 
effort with D. E. Shemansky (regarding the Titan torus data obtained by the 
UVS instrument of the. Voyager spacecraft) was expanded in scope. The first 
step in the expanded collaboration with Shemansky has been to begin to update 
the Titan model of Smyth (1981) to include, the lifetime processes (electron 
impact ionization, charge exchange, and photoionization) for hydrogen opera- 
tive in Saturn’s magnetosphere. Model calculations of the Lyman-a intensity 
in the magnetosphere produced by H atoms lost by Titan may then be properly 
calculated and compared to the UVS measurements of the Voyager spacecraft. 

Such analysis will allow us to determine if gases from the E-ring ic.y satel- 
lites (for which spatial gas envelopes are illustrated for Tethys and Dione in 
Figure 15) or from the planet are relevant in understanding the Voyager UVS 
data. Specific efforts to accomplish these objectives in the third year are 
discussed below. 

3.2 The Role of Titan and Non-Titan Sources 

Efforts in the past year toward understanding the role of Titan (and 
hence non-Titan sources) in the production of Saturn's hydrogen torus have 
centered about (1) gathering information necessary to specify the plasma 
conditions operative at Saturn, (2) calculating the hydrogen lifetime in the 
planetary magnetosphere, and (3) improving the previous Titan hydrogen model 
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of Smyth (1981) by incorporating this spatially varying sink. Pregross in 
those areas was presented as an exhibit at the American Astronomical Society 
16th annual meeting of the Division for Planetary Science (Shemansky, Smith, 
Smyth and Combi, 1984). Plasma information for the planetary magnetosphere 
has been obtained after considerable processing of Voyager PLS data from 
Edward C. Sittler (Goddard Space Flight Center), who has recently become part 
of this larger collaborative effort. The electron and ton densities and 
temperatures along the spacecraft trajectory have been used to generate a two- 
dimensional (longitudinally symmetric) description of the planetary magneto- 
sphere. Further Voyager PLS data have been used to describe solar wind and 
magnetospheric tail plasma conditions. In the magnetosphere, both a cold and 
a hot component of electrons are Included and both light (11 ) and heavy ions 
(0 + and N + ) are considered. 

Lifetime processes for atomic hydrogen in the Saturnian system are 
summarized in Table 3. Using the Voyager 1 outbound plasma data, the lifetime 
of hydrogen in Saturn's equatorial plane is shown in Figure 16 and compared 
with the solar wind charge exchange and photoionization lifetime of H atoms. 
The lifetime value assumed in the earlier Titan torus model calculations of 
Smyth (1981) is also shown. The lifetime decreases fairly rapidly for depar- 
tures normal to the equatorial plane. With the exception of atomic hydrogen 
located inside a radial distance of 10 Saturn radii (Rg), the lifetime based 
upon the Voyager 1 data is larger than that assumed earlier by SmyLh (1981). 
This implies that the nearly longitudinally symmetric distribution of H atoms 
about Saturn (see. Figure 17) calculated earlier by Smyth (1981) will remain 
and, to first order, will be only altered in its radial profile by the 
spatially varying lifetime for H atoms. This alteration is, however, most 
important since it will allow us to divide the Voyager UVS Lyman-a brightness 
data (see Figure 18) into a Titan and non-Titan H source. 

The full two dimensional lifetime of H atoms in the Saturn environment 
has been incorporated in the earlier model of Smyth (1981). Preliminary 
testing of the model was in progress at the close of this project year. This 
work is to be continued under new NASA support and definitive results are 
anticipated in the near future. 
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LOSS P ROCES SES FOR ATOMIC HYDROGEN IN THE SATURNIAN SYSTEM 
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LYMAN ALPHA INTENSITY (R) 



20 0 20 
DISTANCE FROM SATURN (R s ) 


Figure 18. Model vs. UVS Lyman-a Scan of the Hydrogen Torus. 

The Voyager 1 UVS Lyman-a radial scan and a solid line 
for the simple constant density torus model of Broadfoot 
et al. (1981) is shown. The central peak in the data 
points is due to the convolution of Saturn's Lyman-« disk 
by the. UVS slit. 
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IV. COMETARY ATMOSPHERES 


4 . 1 Introdu ction 

Comet ary research at AER represents a broadening and natural extension of 
our previous modeling of extended satellite atmospheres, that was initially 
undertaken for NASA in 19/9 and 1980, Because of the expertise developed at 
AER in modeling extended satellite atmospheres, which are in a collisionless 
regime but remain in a complex gravitation field, we have addressed that area 
of cometary atmospheres where similar conditions are present. The three 
dimensional particle trajectory model developed at AER earlier is now being 
applied to the extended neutral gas clouds observed in comets (e.g. H, C and 
0). In particular, we have constructed a model for the observed spatial dis- 
tribution of cometary hydrogen. The AER three dimensional particle trajectory 
model is ideally suited to modeling cometary 11 atoms which travel under the 
influence of the solar gravitation field and a variable solar radiation 
pressure field owing to resonance scattering of solar Lyman a photons. 

The general state of extended cometary gas cloud models has been one In 
which sets of simplifying assumptions are made to tailor the models to one 
particular aspect of observational data. The two major and somewhat comple- 
mentary modeling efforts to date have been put forth by Keller and his co- 
workers (Keller and Thomas, 1975, Keller and Meier, 1976) and Festou et al. 
(1979). Keller and co-workers have developed models for the extended H cloud 
(>10 6 km) that approximate the source region as a point. These models, 
however, use Maxwe 1 1 -11 o 1 1 xma n n velocity distributions and various approxima- 
tion schemes to calculate atom trajectories under the influence of Lhe solar 
gravitational field and solar radiation pressure. Festou et al. (1979) have 
used the vectorial model, (Festou 1981 a,b) to describe more realistically the 
source region of H In their analysis of the inner coma ^10^ km) Lyman a 
profiles from Comet Kobayashi-Berger-Milon (1975 IX). Their model assumed a 
steady state production ’-ate and neglected (appropriately for these particular 
observations) the effects of the solar gravity field and solar radiation 
pressure. 


4 . 2 The Particle Trajectory Mode l 

Particle trajectory models are not limited to certain spatial regions, as 
in previous modeling efforts. They can be used to calculate correctly the 
three-dimensional atom trajectories in an arbitrarily large spatial domain 
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( > )O fi kin) and at the* same* time to describe properly the spatially extended 

/ p 

source region (10 k4 - l(r km). 

Particle trajectory models have already been suc»*es.*jfully applied to the 

t r 

analysis of neutral radical distributions (iO H to several times ID 4 km). 

Combi and Delsemme (1980) introduced a Monti* Carlo particle trajectory model 
to treat correctly the isotropic ejection of daughter radicals upon photo- 
dissociation of the parent molecules. This model is based therefore on the 
same accurate spatial distribution and kinematics as the vectorial model o f 
Festou (1981a) which was developed at the same time but independently. It 
differs, though, in that it does not involve the lengthy multiple-numerical 
integration, is easily general izable and provides a two-dimensional descrip- 
tion of the coma to Include properly the effect of solar radiation pressure. 
Combi (1980) applied this model with solar radiation pressure to analyze 
sunward and antisunward profiles of CN in Comets Bennett (1970 II) and West 
(1976 VI). 

The cometary hydrogen cloud model, developed at AER, calculates the atom 
trajectories in three dimensions by integrating Newton's equation with a 
fourth order Runge-Kutta routine and includes both the solar gravitational 
field and variable solar radiation pressure, The characteristics adopted for 
the source region (that is, the. starting location and velocity of the H atom 
with respect to the nucleus) are those appropriate for the two step photodis- 
sociation of l^O as described by the Monte Carlo particle trajectory models of 
Combi and Delsemme (1980). The trajectory calculation explicitly Includes the 
orbital elements of the comet, and the absolute time of the observation. The 
placement of an atom on the sky plane is calculated taking into account the 
divergent lines of sight (parallax) appropriate for the relative comet-sun- 
earth geometry. The decay lifetimes of tue parent molecules for H (1^0 and 
OH) are handled in the manner described by Combi and Delsemme (1980), but the 
decay of H, being much larger and having stronger time dependence through the 
changing heliocentric distance, is calculated explicitly in the trajectory 
time-step routine to properly include its spatially changing lifetime. 

The variable radiation pressure acceleration on H atoms as well as the 
instantaneous emission rate of Lyman-a at the observation snap shot are cal- 
culated in the model from the solar disk averaged Lyman-a profile of Lemaire 
et al. (1978). The solar profile is shown in Figure 19. The Lyman-a absorp- 
tion rate by a cometary H atom is a function of the heliocentric velocity 
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(i.e., the position along the profile) as well as the inverse square of the 
heliocentric distance. 

4 . 3 Model Applic a tions 

As a test of the new hydrogen model, we have begun re-analyzing the 
observed I,yman-<x distribution of Comet Kohoutek ( 1973 XII). These data were 
first analyzed by Meier et al. (1976) with their model that was discussed in 
the previous section. They found they could fit the observed isophotes with 
(tu isotropic point source of three Muxwel 1 -Boltzmann distributions having rins 
velocities of 4, 8 and 20 km s~*. This velocity distribution can roughly be 
understood in terms of H prodc.ution from water photodissociation by solar UV 
radiation (Keller and Meier, 1976; Festal et al., 1979). Table 4 shows the 
branching ratios and . roduct velocities for H and OH corresponding to photo- 
dissociation of H 2 O. The major component is in fact at 20 kn s~*. The 
8 km s" 1 component is now clearly understood to result from the photodisso- 
ciation of OH (Jackson, 1980; Schleicher and A' Hearn, 1983; van Dishoeck and 
Dalgano, 1984). The 4 km s~* component had been suggested by Meier et al . to 
result from thermalization of a fraction of the H which was actually produced 
within the collision zone. This low velocity component of H was required for 
smaller heliocentric distances where the gas production rates were higher and 
corresponding decay scale lengths were shorter. This has been demonstrated 
more conclusively by Kitamura et al. (1983). 

Festou et al. (1979) used a velocity distribution as expected from pho- 
todissociation with narrow peaked components at 8 and 20 km s“* which were 
naturally broadened by the spatially extended isotropic dispersion of the 
vectorial model. We have chosen this type of physically motivated picture as 
a starting point for our modeling analysis. Figure 20 shows a contour plot of 
the Lyman-ct emissions from Comet Kohoutek (1973 XII) on January 8, 1974 from 
the original paper by Meier et al. (1976). 

Instead of assuming a simple r” n law for the production rate of H 2 O, in 
our model we have adopted the vaporization curves of Delsemme (1976) as fitted 
to the visual light curves of Comet Kohoutek (see Figure 21). The parent H 2 O 
molecules are assumed to flow isotropically from a point source with a velo- 
city of 0.58 in km s - *, where rj| is in AU (Delsemme, 1982; Whipple, 
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Table 4 

Photoabsorption of Solar UV Radiation by H,0 Vapor 
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0 4 *0.030 
0 8 *0.58Q 
Q 20 * 0.390 


Figure 20. Hydrogen Lyman-a Rocket; Image of Comet Kohoutek by Meier et al. (1976). 

The observed isophotes are shown by the dark contour lines and the model of Meier 
et al. is shown by the light (smooth) lines. The legend gives numerical param- 
eters for their model. 
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Figure i Light Curve of Comet Kohoutek 1973X11. 

The dependence of the gas production rate for Comet Kohoutek on heliocentric 
distance has been set in the particle trajectory model from the light curves 
compiled by Delsemme (1976). The open circles are pre-perihelion points; the 
stars are post-perihelion. The relative scale has been normalized to the pre- 
perihelion production rate at 1 AU. 
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1980). The lifetimes and initial velocities for the various photodissociation 
products are: 


t H 2 0 - 8.2 x 10 4 r H 2 s 

Tqjj * 2.0 x 10 5 rjj 2 s 

Tjj * 2.0 x 10^ r^ 2 s 

v H 2 o = °* 5 8 rjj -0,5 km s-l 

Vqji — 1.2 km s 

20 km s~ l I1 2 0 + hv + OH + H 

V H 8 km s -1 OH + hv + 0 + H 

The modeled Lyman-a image for this case is shown in Figure 22. A 
comparison of the model with the observation shows that the velocity distri- 
bution is still too peaked at 8 and 20 km s - *; the edges of the paraboloid 
envelopes for the two velocity components are apparent in the model. Thus, a 
broader velocity distribution, as would be expected if some thermalization of 
H atoms occurs, is needed. 

In order to test this idea, we have introduced two triangular velocity 
distributions, centered at 8 and 20 km s - -* - with half-widths equal to the 
central velocity. The sum of these two distributions reproduces the three 
Maxwell-Boltzmann distributions of Meier et al. reasonably well. The results 
of this model are shown in Figure 23. The shape and gradient of the intensity 
field match the image reasonably well. 

The effect of thermalization in particle trajectory models for cometary 
atmospheres is currently under study in a separate NASA project at AER and 
will be incorporated into the extended cloud model in the near future. At 
that point, we will see if a true partially-thermalized distribution can 
explain the observed Lyman-a image. For the present, we have developed a 
hydrogen cloud model which can be used to analyze data where collisions are 
not as important (that is, for the right combination of lower gas production 
rate and larger heliocentric distance). The model is an improvement over 
Keller’s earlier work in that it includes (1) a spatially extended source 
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H LYMAN-ALPHA MODEL (RAYLEIGHS) 



DISTANCE FROM THE NUCLEUS (10**6 KM) 


Figure 22. Cometary Hydrogen Particle Trajectory Model: 

Lyman-a Map for Narrow Peaked Velocity Distribution. 

The model run assumed monoenergetic components for hydrogen atoms at 8 and 20 km s”* 
for photodissociation of OH and H£0 respectively. Contour levels are in kiloRay- 
leights and are the same set as in Figure 20. 
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H LYMAN-ALPHA MODEL (RAYLEIGHS) 



DISTANCE FROM THE NUCLEUS (10**6 KM) 


Figure 23. Cometary Hydrogen Particle Trajectory Model: 

Lyman-a Map for Broadened Velocity Distribution. 

The model run assumed triangular velocity distributions centered at 8 and 20 km s“* 
for photodissociation of OH and H2O respectively. These two distributions resemble 

the three Maxwell-Boltzmann distributions used by Meier et al. (1976). Contour 
levels are in kiloRayleighs and are the same set as in Figure 20. 
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region, (2) explicit calculation of atom trajectories in three dimensions, 
ami (3) explicit calculation of the time and velocity dependent lifetime, 
radiation pressure acceleration and excitation rate at each point in the 
cloud. Keller's models computed these quantities only at the position of the 
nucleus or at the Intersection of a line of sight column with the sky plane. 
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